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ABSTRACT 

Remote  sensing  over  long  path  lengths  has  become  of  greater  interest  in  the  terahertz  frequency  region.  Applications 
such  as  pollution  monitoring  and  detection  of  energetic  chemicals  are  of  particular  interest.  Although  there  has  been 
much  attention  to  atmospheric  effects  over  narrow  frequency  windows,  accurate  measurements  across  a  wide  spectrum 
is  lacking.  The  water  vapor  continuum  absorption  spectrum  was  investigated  using  Fourier  Transform  Spectroscopy. 

The  continuum  effect  gives  rise  to  an  excess  absorption  that  is  unaccounted  for  in  just  a  resonant  line  spectrum 
simulation.  The  transmission  of  broadband  terahertz  radiation  from  0.300THz  -  1.5THz  through  air  with  varying  relative 
humidity  levels  was  recorded  for  multiple  path  lengths.  From  these  data,  the  absorption  coefficient  as  a  function  of 
frequency  was  determined  and  compared  with  model  calculations.  The  intensity  and  location  of  the  strong  absorption 
lines  were  in  good  agreement  with  spectral  databases  such  as  the  2008  FIITRAN  database  and  the  JPL  database. 
However,  a  noticeable  continuum  effect  was  observed  particularly  in  the  atmospheric  transmission  windows.  A  small 
discrepancy  still  remained  even  after  accounting  for  continuum  absorption  using  the  best  available  data  from  the 
literature.  This  discrepancy,  when  projected  over  a  one  kilometer  path  length,  typical  of  distances  used  in  remote 
sensing,  can  cause  a  30dB  difference  between  calculated  and  observed  attenuation.  From  the  experimental  and  resonant 
line  simulation  spectra  the  air-broadening  continuum  parameter  was  calculated  and  compared  with  values  available  in 
the  literature. 
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1  INTRODUCTION 

The  propagation  of  electromagnetic  radiation,  in  particular  at  terahertz  frequencies,  through  the  Earth’s  atmosphere  is  an 
important  field  of  study.  Applications  like  remote  sensing,  imaging,  detection  of  concealed  items,  astronomical 
observations,  and  more  are  used  in  industries  like  defense,  medical,  environmental,  and  academia  and  are  more 
commonly  employing  terahertz  frequencies.  Water  vapor  has  thousands  of  rotational  and  vibrational  absorption  lines 
positioned  throughout  the  whole  electromagnetic  spectrum  and  is  a  main  contributor  to  atmospheric  attenuation  in  the 
terahertz  regime.  Propagation  of  radiation  over  extended  distances  is  not  feasible  at  many  frequencies.  Measurements  to 
be  performed  in  this  frequency  region  require  significant  attention  to  the  specific  environment. 

The  absorption  peaks  of  water  vapor  have  been  extensively  studied1'3  however  the  window  regions  between  the  strong 
absorption  lines  are  more  difficult  to  characterize.  Current  models4'8  require  the  use  of  an  empirical  continuum 
contribution  to  the  spectrum  for  an  accurate  comparison  to  experimental  spectra.  In  the  microwave  region  there  is  much 
evidence  that  the  continuum  contribution  scales  as  the  square  of  the  frequency  with  a  negative  dependence  on 
temperature.8"10  However,  at  higher  frequencies  in  the  terahertz  region  only  a  few  studies  have  seen  the  quadratic 
dependence  on  frequency.11"13  Further  complicating  the  issue,  there  is  uncertainty  as  to  the  cause  of  the  continuum  as 
well  as  the  values  of  the  continuum  parameters.  Since  the  continuum  contribution  is  empirical,  it  depends  on  both 
experimental  data  as  well  as  model  predictions.  The  choice  of  line  parameters,  line  shape  function,  and  number  of 
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absorption  lines  used  in  the  model  can  significantly  affect  the  calculation  of  the  continuum  parameters.  Also  a  limited 
number  of  experimental  studies  performed  at  terahertz  frequencies  restricts  a  comprehensive  and  accurate  determination 
of  the  continuum  parameters.  Only  a  few  broadband  studies  have  been  performed  in  the  terahertz  region,  most  of  which 
were  performed  in  a  laboratory  environment  using  either  pure  nitrogen  or  oxygen  gas  as  the  foreign  broadening  gas.1215 
Control  over  the  experimental  conditions  is  gained  through  the  sacrifice  of  the  broadening  gas,  which  is  of  a  different 
composition  from  Earth’s  atmosphere.  Other  works  have  studied  atmospheric  air,1 1  1618  however  these  studies  do  not 
have  control  over  many  of  the  experimental  parameters. 

The  present  study  used  dry  air  as  the  foreign  broadening  gas,  which  has  the  same  composition  of  atmospheric  air,  while 
maintaining  control  over  the  experimental  parameters.  Broadband  transmission  data  was  taken  between  0.3-1.5THz  over 
multiple  path  lengths  and  humidity  levels  to  determine  the  absorption  coefficient  of  atmospheric  air.  These  results  were 
combined  with  simulated  absorption  coefficient  data  calculated  using  line  parameters  from  the  2008  HITRAN 
Database19  (HITRAN)  to  find  the  air-broadened  continuum  coefficient.  When  the  continuum  contribution  is 
incorporated  into  the  absorption  coefficient  calculations,  the  resulting  spectrum  can  be  used  to  identify  the  best 
frequency  range  for  which  an  experiment  is  to  be  performed. 


2  PROCEDURE 


2.1  Experimental  setup 

Fourier  Transform  Spectroscopy  (FTS)  was  used  to  collect  transmission  data  of  terahertz  radiation  through  air  with 
varying  relative  humidity  levels  and  path  lengths.  A  Pike  Technologies  variable  path  length  gas  cell  was  used  to  contain 
the  sample  and  allow  for  multiple  path  length  measurements.  Figure  1  shows  a  schematic  of  the  experimental  setup 
along  with  a  diagram  of  the  variable  path  length  gas  cell.  The  cell  was  mounted  into  the  sample  chamber  of  a  Bruker 
Vertex  80V  Fourier  Transform  Infrared  (FTIR)  Spectrometer  with  a  mercury  arc  lamp  source  aligned  to  the  optical  path 
of  the  spectrometer.  The  transmitted  radiation  was  collected  using  a  liquid  helium  cooled  bolometer  from  IR  labs  with  a 
vacuum  pump  attached  to  the  liquid  helium  chamber  to  lower  the  operating  temperature  to  1.6K. 


The  variable  path  length  cell  used  in  the  experiment  was  able  to  achieve  l-6m  path  lengths  in  lm  increments.  The  path 
length  was  adjusted  through  a  micrometer  screw  on  the  top  of  the  cell  without  the  need  to  realign  the  setup.  The  cell  was 
connected  to  a  vacuum  pump  and  sample  containers.  The  pressure  of  the  system  was  monitored  using  external 
capacitance  manometers  allowing  for  measurements  from  O.Ol-lOOOTorr.  In  addition  a  Honeywell  HIH-4000  humidity 
sensor  was  mounted  in  the  cell,  which  allowed  for  an  independent  measure  of  the  relative  humidity,  and  hence  the  partial 
pressure  of  water  vapor,  inside  the  cell. 
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2.2  Transmission  measurements  and  absorption  coefficient  calculations 

Transmission  data  were  recorded  from  300-1500GHz  with  a  resolution  of  3GHz.  The  FTIR  used  a  125pm  Mylar 
beamsplitter,  which  has  a  null  at  720GHz.  The  bolometer  also  had  a  cold  filter,  which  attenuated  radiation  above 
1500GHz.  As  stated  previously,  each  dataset  had  6  transmission  scans,  one  for  each  path  length.  A  total  of  4  relative 
humidity  levels  were  studied,  17.89,  36.79,  48.03,  and  70.84%  relative  humidity.  These  relative  humidity  levels 
correspond  to  3.60,  8.04,  10.05,  and  15.23Torr  partial  pressure  of  water  vapor  respectively  at  296K. 

To  collect  a  transmission  scan  the  cell  was  evacuated  to  below  lOmTorr  and  sealed  for  background  measurements. 
Background  measurements  were  recorded  for  l-6m  path  lengths  consecutively.  After  the  background  scans  were 
completed,  distilled  water  vapor  followed  by  dry  air  was  introduced  into  the  cell  and  then  sealed.  Sample  scans  for  path 
lengths  of  1  -6m  were  then  completed  and  normalized  to  the  corresponding  background  scan. 

The  absorption  coefficient  was  then  calculated  from  the  transmission  data  using  the  Beer-Lambert  Law,  x=exp(-al), 
where  x  is  the  normalized  transmission,  a  is  the  absorption  coefficient,  and  1  is  the  path  length  of  the  cell.  The 
absorption  coefficient  is  then  found  to  be  the  slope  of  the  straight  line  fit  of  the  logarithm  of  the  transmission  versus  the 
path  length,  as  shown  in  Eq.  (  1  ). 

al  =  -ln(r)  ( 1 ) 


3  ABSORPTION  COEFFICIENT  CAECUEATIONS 

A  detailed  description  of  the  absorption  coefficient  calculations  can  be  found  elsewhere,20  so  only  a  brief  description  of 
the  calculations  will  be  given  here.  The  calculations  for  an  absorption  coefficient  spectrum  were  performed  in  five 
successive  steps  for  direct  comparison  with  experimental  datasets. 

Step  one  of  the  calculations  is  to  simulate  the  self-broadened  absorption  line.  The  line  parameters  for  each  absorption 
line  were  obtained  from  HITRAN  and  were  used  in  Eq.(  2  ),the  Van  Vleck-Weisskopf  line  shape  function.21 


(2) 


In  Eq.  (  2  ),  and  a^max  are  the  absorption  coefficient  of  the  kth  line  in  the  database  and  the  amplitude  of  the  self 
broadened  absorption  line  in  units  of  m_1  and  v,  Vk,  and  Avk  are  the  frequency  of  the  incident  radiation,  line  center 
frequency,  and  line  half  width  respectively  in  MHz. 

Step  two  of  the  calculations  is  to  account  for  the  effects  of  foreign  gas  broadening.  Equation  (  2  )  is  again  used  to 
calculate  the  foreign  gas  broadened  line  shape.  The  whole  line  shape  is  then  multiplied  by  a  normalization  factor  so  that 
the  integrated  line  strength  is  the  same  as  the  self  broadened  line  shape  from  step  one.  Figure  2  shows  a  comparison  of 
the  two  pressure-broadened  line  shapes  as  well  as  a  Doppler-broadened  line  shape.  As  can  be  seen  in  the  figure,  the 
effect  of  foreign  gas  broadening  is  to  decrease  the  peak  amplitude  and  increase  the  line  width. 

Step  three  is  to  take  the  sum  of  all  of  the  foreign  gas  broadened  absorption  lines  in  the  frequency  region  to  produce  the 
resonant  line  absorption  coefficient  spectrum.  Step  four  is  then  to  account  for  the  continuum  contribution  to  the 
absorption  coefficient.  The  water  vapor  continuum  can  be  accurately  modeled  by  Eq.  (  3  )  "10,13"14 

«,  M  -  v2P,  (p,c,8"‘  +  Pffi"’ )  (3) 
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Frequency  (GHz) 


Figure  2.  A  Doppler-broadened  absorption  line  with  the  partial  pressure  of  water  vapor  equal  to  the  critical  pressure  (black)  and  two 
pressure-broadened  absorption  lines;  a  self-broadened  line  with  15.23Torr  of  water  vapor  (red)  and  a  foreign  gas  broadened  line  with 
15.23Torr  of  water  vapor  and  746.77Torr  of  air  (blue). 


In  Eq.  ( 3  ),  ac  is  the  continuum  contribution  to  the  absorption  coefficient  in  units  of  m'1,  Cs  and  Cf  are  the  self  and 
foreign  continuum  coefficients  in  units  of  m"V(Torr2MHz2),  Ps  and  Pf  are  the  sample  and  foreign  gas  pressures  in  units  of 
Torr,  the  exponential  terms  ns  and  nf  are  the  self  and  foreign  temperature  exponents,  and  0  is  a  reference  temperature 
divided  by  the  gas  temperature  in  Kelvin. 

The  continuum  component  of  the  absorption  coefficient,  as  calculated  from  Eq.  (  3  ),  is  then  added  to  the  resonant  line 
absorption  coefficient  spectrum  calculated  in  step  three  to  produce  the  total  absorption  coefficient.  Step  five  is  to 
account  for  the  instrumentation  distortion  of  the  experimental  apparatus.  The  FTIR  uses  discrete  resolution  bins  of  finite 
width  in  the  Fast  Fourier  Transform  of  the  transmission  spectra.  This  distortion  is  modeled  by  applying  the  same 
discrete  resolution  bins  to  the  calculated  spectra. 

The  effects  of  steps  3-5  and  the  resulting  spectra  can  be  seen  in  Figure  3.  The  continuum  absorption  has  a  small  but 
measurable  contribution  to  the  total  absorption  coefficient.  This  contribution,  however,  is  only  significant  in  the 
transmission  windows  between  the  strong  spectral  lines.  At  the  peaks  of  the  strong  spectral  lines  the  continuum 
contribution  to  the  absorption  coefficient  is  unnoticeable.  In  the  figure  this  can  be  clearly  seen  as  a  difference  between 
the  green  and  red  spectra  in  the  windows,  while  at  the  peaks  the  green  and  red  spectra  overlap  each  other  and  are 
indistinguishable.  Also  shown  in  the  Figure  3  are  the  effects  of  the  instrumentation  distortion.  The  effect  is  to  lower  and 
broaden  the  spectral  lines.  The  distortion  is  only  significant  in  the  strong  spectral  lines  and  is  unnoticeable  in  the 
transmission  windows  because  the  absorption  coefficient  has  a  much  smaller  gradient  in  these  regions.  Again  this  can  be 
clearly  seen  in  the  figure,  however  for  the  distortion  the  red  and  blue  spectra  overlap  and  are  indistinguishable  in  the 
window  regions  and  in  the  peaks  there  is  a  clear  difference  between  the  red  and  blue  spectra. 

It  is  necessary  to  account  for  any  instrumentation  distortion  in  the  calculated  spectrum  to  make  an  accurate  comparison 
with  the  experimental  peak  amplitudes.  The  difference  between  the  peak  amplitudes  of  a  calculated  spectrum  with  and 
without  instrumentation  distortion  for  a  resolution  of  3GHz  can  be  significant.  If  this  distortion  is  unaccounted  for  direct 
comparison  of  calculated  spectra  and  experimental  spectra  will  yield  poor  results.  This  distortion  also  restricts  the 
comparison  of  two  experimental  spectra  directly  without  similar  instrumentation  distortion  or  adjustments  to  one  or  both 
of  the  spectra. 
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Figure  3.  The  results  from  the  absorption  coefficient  calculations.  The  air-broadened  spectrum  (green),  continuum  absorption  (pink), 
air-broadened  spectrum  with  continuum  absorption  (dashed  red),  and  air-broadened  spectrum  with  continuum  absorption  and 
instrumentation  distortion  (blue).  All  plots  are  15.23Torr  of  water  vapor  and  746.77Torr  or  air. 


4  RESULTS 

The  transmission  plots  of  l-6m  path  lengths  for  terahertz  radiation  through  air  with  relative  humidity  levels  of  17.89, 
36.78,  48.03,  and  70.84%  were  recorded  and  the  dataset  for  70.84%  relative  humidity  is  shown  in  Figure  4.  The  width  of 
the  studied  frequency  region  was  limited  by  the  power  emitted  from  the  source  and  a  cold  filter  in  the  bolometer  on  the 
low  and  high  end  of  the  region.  As  seen  in  the  figure,  there  is  a  region  of  significant  noise  centered  at  720GHz,  which  is 
due  to  a  beamsplitter  null.  Also  seen  in  the  figure  are  saturation  effects  in  many  of  the  strong  absorption  lines, 
particularly  in  the  frequency  range  from  1.05-1.25THz. 

The  absorption  coefficient  was  calculated  for  each  dataset  according  to  the  process  outline  in  Section  2.2  and  is  shown 
for  70.84%  relative  humidity  in  Figure  5.  Slocum  et  al.20  published  the  full  results  of  the  study  including  plots  of  all  four 
datasets  for  the  transmission  and  absorption  coefficient.  The  effects  of  saturation  can  be  seen  in  the  strong  absorption 
lines  as  regions  of  missing  data.  Saturation  limited  the  calculation  of  the  absorption  coefficient  to  values  less  than  4m"1. 
The  statistical  error  for  each  data  point  was  calculated  from  the  fitting  routine  of  the  absorption  coefficient.  The  average 
error  across  the  data  range,  excluding  regions  of  low  signal  to  noise,  was  8.85,  4.00,  3.67,  and  2.87%  for  the  17.89, 

36.78,  48.03,  and  70.84%  relative  humidity  datasets  respectively. 
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Figure  4.  The  experimental  transmission  of  air  with  a  relative  humidity  of  70.84%  for  l-6m.  The  red,  orange,  yellow,  green,  blue, 
and  purple  plots  correspond  to  l-6m  respectively. 


The  empirical  continuum  absorption  data  was  found  from  taking  the  difference  of  the  experimental  absorption 
coefficient  data  and  the  calculated  absorption  coefficient  data  without  compensating  for  continuum  absorption.  This  was 
then  used  to  determine  the  air-broadened  continuum  coefficient  by  fitting  the  continuum  absorption  to  Eq.  (  3  ).  The 
fitting  parameter  was  Cf*=Cf0”f  since  all  of  the  data  from  the  present  study  were  at  a  constant  temperature  of  296K  so  0nf 
is  a  multiplicative  constant  for  the  present  study.  The  value  obtained  of  Cf*=1.158X10"18m"V(Torr2MHz2)  was  found  to 
minimize  the  mean  square  error  of  the  four  datasets.  The  fitting  routine  employed  the  value  of  Cs*  from  Podobedov  et 
al. 13-14  of  1.76X10"17m"V(Torr2MHz2),  which  yielded  a  minimum  for  the  mean  square  error  when  compared  with  other 
values  from  the  literature. 


Figure  5.  The  experimentally  determined  absorption  coefficient  for  air  with  a  relative  humidity  of  70.84%. 
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5  ANALYSIS 

The  calculated  and  experimental  absorption  coefficients  were  converted  from  units  of  m'1  to  dB/km  in  the  remainder  of 
the  paper  for  ease  of  use*.  The  experimental  spectra  were  compared  with  calculated  values  from  the  process  outlined  in 
Section  3.  The  continuum  contribution  to  the  absorption  coefficient  was  calculated  using  parameter  sets  from  Kuhn  et 
al.,9  Podobedov  et  al.,1314  and  the  present  study.  The  residual  of  each  of  these  calculated  spectra  with  the  experimental 
spectra  was  a  minimum  for  the  continuum  parameter  set  from  the  present  study.  The  results  of  this  comparison  can  be 
seen  in  Figure  6.  Given  these  results,  only  the  continuum  parameter  set  from  the  present  study  is  used  for  comparison  to 
the  experimental  data.  Figure  7  shows  the  calculated  spectrum  along  with  the  experimental  absorption  coefficient 
spectrum.  As  can  be  seen,  when  the  continuum  is  unaccounted  for  there  are  significant  discrepancies  between  the 
experimental  and  calculated  spectra  in  the  window  regions.  However,  after  the  continuum  contribution  has  been  added 
there  is  very  good  agreement  between  the  experimental  and  calculated  spectra.  It  can  also  be  seen,  as  stated  previously, 
that  in  the  peaks  of  the  spectra  there  is  no  noticeable  difference  between  the  two  calculated  spectra  and  the  continuum 
absorption  only  makes  a  significant  contribution  in  the  window  regions. 


Figure  6.  A  histogram  of  the  percent  difference  with  a  bin  size  of  10  points  between  the  observed  spectra  and  the  calculated  spectra 
using  different  continuum  parameter  sets.  The  plot  shows  comparisons  with  calculations  using  the  continuum  parameter  set  from  the 
present  study  (red),  Kuhn  et  al.  (green),  Podobedov  et  al.  (blue)  and  with  no  continuum  (black). 


Figure  8  shows  an  expanded  view  of  the  atmospheric  window  at  870GHz.  As  discussed  previously,  the  calculated 
spectrum  without  continuum  effects  does  a  poor  job  of  reproducing  the  experimental  spectrum.  However,  when 
continuum  absorption  is  accounted  for  in  the  calculations  there  is  good  agreement  between  the  calculated  and 
experimental  spectra.  Also  seen  in  the  figure  is  an  absorption  peak  located  at  916GHz.  It  can  be  seen  that  for  the 
absorption  line  it  is  not  necessary  to  account  for  continuum  absorption  to  yield  and  accurate  comparison,  both  calculated 
spectra  are  in  good  agreement  wit  the  experimental  spectrum  at  the  absorption  line. 


*  lm_1=4. 343X1 03dB/km 
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Figure  7.  A  comparison  of  the  experimental  absorption  coefficient  (red)  and  the  calculated  absorption  coefficient  both  with  (green) 
and  without  (black)  a  continuum  contribution.  All  plots  are  for  10.05Torr  of  water  vapor  and  750.95Torr  of  air 
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The  calculation  process  outlined  in  Section  3  was  also  used  to  model  a  broad  range  of  experimental  absorption  data 
found  in  the  literature,  the  results  of  which  can  be  seen  in  Figure  9.  The  calculations  provide  good  agreement  with  the 
experimental  data  collected  using  four  different  collection  techniques;  cavity  Q  measurements,  ’  ’  Time  Domain 
Spectroscopy,26  direct  attenuation  measurements20  and  FTS  measurements  from  the  present  study.  The  calculations  are 
also  able  to  accurately  model  a  wide  range  of  conditions  including  self-broadening  at  low  pressures,  and  foreign  gas 
broadening  at  both  low  and  high  water  vapor  pressures  for  temperatures  ranging  from  296-376K.  The  calculations  for 
this  figure  used  the  continuum  parameter  set  from  Kuhn  et  al.  The  experimental  data  displayed  in  the  figure  were  taken 
at  varying  temperatures.  The  present  study  was  not  able  to  determine  the  temperature  dependence  of  the  continuum, 
which  is  necessary  to  accurately  model  the  wide  range  of  data  in  the  figure.  Another  important  consideration  when 
comparing  experimental  data  obtained  using  different  measurement  techniques  is  that  each  technique  will  have  a 
different  instrumentation  distortion  on  the  spectra.  Therefore  instrumentation  distortion  must  be  taken  into  account 
when  comparing  two  experimental  datasets  or  an  experimental  with  a  calculated  dataset.  These  considerations  were 
taken  into  account  in  Figure  9.  When  applied  correctly  these  calculations  can  be  very  useful  in  setting  up  an  experiment 
or  designing  a  device  that  will  incorporate  the  propagation  of  terahertz  radiation  through  the  open  atmosphere.  Given 
some  range  of  operating  conditions  the  calculations  can  establish  how  much  signal  can  be  expected  at  the  receiver  after 
travelling  some  know  distance.  This  knowledge  can  aid  in  the  choice  of  the  source  and  receiver  used  as  well  as  setting 
realistic  propagation  distances. 


Calculations  (dB/km) 

Figure  9.  Correlation  Plot  of  experimental  data  from  the  literature  and  the  absorption  coefficient  calculations. 


Figure  10  shows  a  comparison  of  calculated  absorption  coefficient  values  for  ITorr  of  water  vapor  and  750Torr  of 
nitrogen  for  296K  at  two  different  frequencies.  The  spectra  have  a  continuum  absorption  contribution  and  have  the  same 
calculation  parameters  except  for  the  high  frequency  limit.  As  can  be  seen  in  the  figure,  the  absorption  coefficient 
increase  with  the  increase  of  the  high  frequency  limit.  This  is  expected  as  a  larger  high  frequency  limit  means  the 
inclusion  of  more  absorption  lines  each  contributing  to  the  absorption  coefficient  at  every  frequency.  However,  the 
effect  reaches  saturation  when  the  lines  are  sufficiently  far  from  the  studied  frequency.  The  frequency  difference 
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required  to  be  sufficiently  far  may  not  be  a  constant  or  even  a  proportional  value,  but  is  influenced  by  the  magnitude  of 
the  absorption  coefficient  at  the  studied  frequency.  The  percent  change  in  the  absorption  coefficient  for  a  given  change 
in  the  high  frequency  limit  is  larger  when  the  magnitude  of  the  absorption  coefficient  is  larger.  The  reference  value  for 
the  846GHz  data  is  1.8977X10"3m_l  and  for  the  1.56THz  data  it  is  1.497 1X1  O^m"1.  The  percent  change  in  the  absorption 
coefficient  is  smaller  for  the  larger  of  the  two  reference  values.  This  is  because  the  additional  absorption  provided  by  the 
additional  lines  is  on  the  same  order  of  magnitude  for  both  frequencies.  Therefore  a  smaller  reference  value  will  yield  a 
larger  percent  change.  Since  frequencies  below  ITHz  tend  to  have  smaller  absorption  coefficient  values,  the  choice  of 
the  high  frequency  limit  influences  the  calculations  more  at  these  frequencies.  A  calculation  to  3.56THz  is  sufficient  to 
reach  the  saturation  effect  at  1.56THz,  while  for  a  measurement  at  846GHz  a  calculation  to  4.46THz  is  required  to  reach 
saturation. 


A  846GHz 
•  1.560THZ 

▲ 


10  - 


I 

o 


8  - 


6  - 


2  - 


1  2  3  4  5  6  7  8 

Frequency  beyond  data  point  included  in  calculation  (THz) 

Figure  10.  Comparison  of  calculated  absorption  coefficients  with  different  high  frequency  limits.  The  reference  absorption 
coefficient  value  for  the  percent  difference  calculation  is  from  the  calculation  with  the  high  frequency  limit  750GHz  beyond  the 
specified  frequency. 


6  ABSORPTION  COEFFICIENT  PREDICTIONS 

The  calculation  routine  described  in  Section  3  was  used  to  create  predictions  of  the  absorption  coefficient  of  atmospheric 
air  for  varying  relative  humidity  levels.  Datasets  were  created  for  0-3THz  at  296K  with  10-90%  relative  humidity  in 
10%  increments.  The  calculations  contain  a  continuum  absorption  contribution  calculated  using  the  parameter  set  from 
the  present  study  and  do  not  account  for  any  instrumentation  distortion.  Figure  1 1  shows  four  such  plots  for  10,  30,  50 
and  90%  relative  humidity.  Marked  in  the  figure  are  16  atmospheric  windows  where  the  absorption  coefficient  has  been 
determined  and  displayed  in  Table  1. 
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Frequency  (GHz) 

Figure  11.  Calculated  absorption  coefficient  values  of  water  vapor  at  296K  for  10,  30,  50,  and  90  percent  relative  humidity. 
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Table  1 .  A  list  of  the  calculated  absorption  coefficient  values  in  dB/km  of  atmospheric  air  of  varying  relative  humidity  at  296K  for 
the  windows  identified  in  Figure  11. 


Frequency  (GHz) 

10% 

20% 

30% 

40% 

50% 

60% 

70% 

80% 

90% 

215.543 

1.170 

2.370 

3.630 

4.940 

6.290 

7.690 

9.140 

10.60 

12.20 

340.630 

3.250 

6.610 

10.10 

13.80 

17.50 

21.40 

25.40 

29.50 

33.80 

408.657 

5.520 

11.20 

17.20 

23.30 

29.60 

36.10 

42.90 

49.80 

57.00 

668.519 

15.30 

31.00 

47.40 

64.30 

81.70 

99.50 

118.0 

137.0 

156.0 

845.159 

18.20 

37.10 

56.80 

77.20 

98.30 

120.0 

143.0 

166.0 

190.0 

1023.632 

71.30 

144.0 

220.0 

297.0 

377.0 

458.0 

542.0 

628.0 

716.0 

1348.493 

79.30 

161.0 

245.0 

333.0 

422.0 

514.0 

609.0 

707.0 

807.0 

1501.388 

76.90 

156.0 

238.0 

324.0 

411.0 

502.0 

595.0 

691.0 

790.0 

1829.587 

252.0 

511.0 

778.0 

1050 

1330 

1620 

1920 

2220 

2540 

1979.331 

140.0 

285.0 

435.0 

589.0 

749.0 

913.0 

1080 

1260 

1440 

2105.282 

230.0 

467.0 

711.0 

963.0 

1220 

1490 

1760 

2040 

2330 

2304.914 

385.0 

781.0 

1190 

1610 

2040 

2480 

2940 

3410 

3890 

2522.282 

253.0 

513.0 

784.0 

1060 

1350 

1650 

1960 

2270 

2600 

2715.789 

648.0 

1310 

1990 

2700 

3420 

4150 

4910 

5690 

6480 

2839.913 

515.0 

1040 

1590 

2150 

2730 

3320 

3930 

4550 

5190 

2916.306 

615.0 

1250 

1900 

2570 

3260 

3970 

4690 

5440 

6200 

7.  Conclusions 

Presented  in  this  paper  are  broadband  absorption  measurements  of  air  with  varying  levels  of  relative  humidity  for  the 
0. 300-1. 500THz  frequency  regime.  FTS  measurements  of  transmission  spectra  were  recorded  over  multiple  path  lengths 
and  used  to  determine  the  absorption  coefficient  as  a  function  of  frequency.  These  experimental  absorption  spectra 
along  with  calculated  absorption  spectra  were  used  to  calculate  the  air-broadened  continuum  coefficient  to  be 
1.158X10'lsm"7(Torr2MHz2).  This  value  is  within  the  scatter  of  values  previously  reported  in  the  literature.  The 
calculation  routine  was  also  used  to  model  experimental  water  vapor  absorption  data  from  other  authors.  The 
calculations  provide  good  agreement  to  the  data  for  four  measurement  techniques  with  a  wide  range  of  experimental 
parameters  over  8  orders  of  magnitude  in  the  absorption  coefficient.  This  calculation  routine  can  be  used  as  an  accurate 
predictive  tool  for  atmospheric  absorption  of  terahertz  radiation  when  the  continuum  contribution  is  accounted  for.  It 
was  shown  that  instrumentation  distortion  must  be  accounted  for  before  a  comparison  of  data  can  be  made  between 
different  collection  techniques.  If  the  continuum  is  not  included,  the  discrepancy  between  the  calculations  and 
experiment  could  be  as  large  as  90dB  over  a  one  kilometer  path  length. 
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